High spatial entanglement via chirped quasi-phase-matched optical parametric 

down-conversion 
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By making use of the spatial shape of paired photons, parametric down-conversion allows the gen- 
eration of two-photon entanglement in a multidimensional Hilbert space. How much entanglement 
can be generated in this way? In principle, the infinite-dimensional nature of the spatial degree 
of freedom renders unbounded the amount of entanglement available. However, in practice, the 
specific configuration used, namely its geometry, the length of the nonlinear crystal and the size 
of the pump beam, can severely limit the value that could be achieved. Here we show that the 
use of quasi-phase-matching engineering allows to increase the amount of entanglement generated, 
reaching values of tens of ebits of entropy of entanglement under different conditions. Our work 
thus opens a way to fulfill the promise of generating massive spatial entanglement under a diverse 
variety of circumstances, some more favorable for its experimental implementation. 

PACS numbers: 03.67.Bg, 03.65. Aa, 42.50.Dv, 42.65.Lm 



Entanglement is a genuine quantum correlation be- 
tween two or more parties, with no analogue in classi- 
cal physics. During last decades it has been recognized 
as a fundamental tool in several quantum information 
protocols, such as quantum teleportation [l|, quantum 
cryptography Q and quantum key distribution and 
distributed quantum computing [J]. 

Nowadays, spontaneous parametric down-conversion 
(SPDC), a process where the interaction of a strong 
pump beam with a nonlinear crystal mediates the emis- 
sion of two lower-frequency photons (signal and idler), is 
a very convenient way to generate photonic entanglement 
Photons generated in SPDC can exhibit entangle- 
ment in the polarization degree of freedom Q , frequency 
and spatial shape @, Q . One can also make use of a 
combination of several degrees of freedom 1(| U | ■ 

Two-photon entanglement in the polarization degree of 
freedom is undoubtedly the most common type of gener- 
ated entanglement, due both to its simplicity, and that 
it suffices to demonstrate a myriad of important quan- 
tum information applications. But the amount of en- 
tanglement is restricted to 1 ebit of entropy of entangle- 
ment [l2| , since each photon of the pair can be generally 
described by the superposition of two orthogonal polar- 
izations (two-dimensional Hilbert space). On the other 
hand, frequency and spatial entanglement occurs in an 
infinite dimensional Hilbert space, offering thus the pos- 
sibility to implement entanglement that inherently lives 
in a higher dimensional Hilbert space (qudits). 

Entangling systems in higher dimensional systems (fre- 
quency and spatial degrees of freedom) is important both 
for fundamental and applied reasons. For example, noise 
and decoherence tend to degrade quickly quantum corre- 



lations. However, theoretical investigations predict that 
physical systems with increasing dimensions can main- 
tain non-classical correlations in the presence of more 
hostile noise 13|, ll4j . Higher dimensional states can also 
exhibit unique outstanding features. The potential of 
higher-dimensional quantum systems for practical appli- 
cations is clearly illustrated in the demonstration of the 
so-called quantum coin tossing, where the power of higher 
dimensional spaces is clearly visible [IH • 

The amount of spatial entanglement generated de- 
pends of the SPDC geometry used (collinear vs non- 
collinear), the length of the nonlinear crystal (L) and 
the size of the pump beam (wq). To obtain an initial es- 
timate, let us consider a collinear SPDC geometry. Un- 
der certain approximations jl6|, the entropy of entan- 
glement can be calculated analytically. Its value can be 
shown to depend on the ratio L/Ld, where — k p ur^/2 
is the Rayleigh range of the pump beam and k p is its 
longitudinal wavenumbcr. Therefore, large values of the 
pump beam waist wq and short crystals are ingredients 
for generating high entanglement [17] . However, the use 
of shorter crystals also reduce the total flux-rate of gen- 
erated entangled photon pairs. Moreover, certain ap- 
plications might benefit from the use of focused pump 
beams. For instance, for L = 1 mm, wg = 200/xm and 
kp = 15.7/im -1 , one obtains E ~ 9 [16]. For a longer 
crystal of L = 20 mm, the amount of entanglement is 
severely reduced to E ~ 5 ebits. 

We put forward here a scheme to generate massive spa- 
tial entanglement, i. e., an staggering large value of the 
entropy of entanglement, independently of some relevant 
experimental parameters such as the crystal length or 
the pump beam waist. This would allow to reach even 
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FIG. 1. Scheme of SPDC in a linearly chirped quasi-phase- 
matched nonlinear crystal. The pump beam is a Gaussian 
beam, and p and q designate the transverse wave numbers of 
the signal and idler photons, respectively. Ko is the grating 
wave-vector at the input face of the nonlinear crystal, and 
Ko — Q.L at its output face. The signal and idler photons can 
have different polarizations or frequencies. The different col- 
ors (or different direction of arrows) represent domains with 
different sign of the nonlinear coefficient. 



larger amounts of entanglement that possible nowadays 
with the usual configurations used, or to attain the same 
amount of entanglement but with other values of the 
nonlinear crystal length or the pump beam waist better 
suited for specific experiments. 

Our approach is based on a scheme originally used to 
increase the bandwidth of parametric down-conversion 
18 20]. A schematic view of the SPDC configuration 
is shown in FigJT] It makes use of chirped quasi-phase- 
matching (QPM) gratings with a linearly varying spatial 
frequency given by K g (z) — Kq — a(z + L/2), where Kq 
is the grating's spatial frequency at its entrance face (z = 
—L/2), and a is a parameter that represents the degree 
of linear chirp. The period of the grating at distance z is 
p(z) = 2n/ Kg(z), so that the parameter a writes 



2tt p f - pi 
L PfPi 



(1) 



where Pi is the period at the entrance face of the crystal, 
and pf at its output face. 

The key idea is that at different points along the non- 
linear crystal, signal and idler photons with different fre- 
quencies and transverse wavenumbers can be generated, 
since the continuous change of the period of the QPM 
gratings allows the fulfillment of the phase-matching con- 
ditions for different frequencies and transverse wavenum- 
bers. If appropriately designed narrow-band interference 
filters allow to neglect the frequency degree of freedom 
of the two-photon state, the linearly chirped QPM grat- 
ing enhance only the number of spatial modes generated, 
leading to a corresponding enhancement of the amount 
of generated spatial entanglement. 





FIG. 2. Weight of the Schmidt coefficients X„i for (a) a = 
fj,rn~ 2 and (b) a = 10 x 10~ 6 fj,rn~ 2 . The nonlinear crystal 
length is L = 20 mm and the pump beam waist is wo = 
200 urn. 



In order to determine how much spatial entanglement 
can be generated in SPDC with the use of chirped QPM, 
let us consider a nonlinear optical crystal illuminated 
by a quasi-monochromatic laser Gaussian pump beam 
of waist wq. Under conditions of collinear propagation 
of the pump, signal and idler photons with no Poynting 
vector walk-off, which would be the case of a noncrit- 
ical type-II quasi-phase matched configuration, the am- 
plitude of the quantum state of the generated two-photon 
pair generated in SPDC reads in transverse wavenumber 
space 



I*) 



dpdq$ (p,q) |p) s |q)i, 



(2) 
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where p (q) is the transverse wavenumber of the sig- 
nal (idler) photon. \P is the joint amplitude of the two- 
photon state, so that (p, q) | 2 is the probability to de- 
tect a signal photon with transverse wave-number p in 
coincidence with an idler photons with q. 



The joint amplitude that describes the quantum state 
of the paired photons generated in a linearly chirped 
QPM crystal, using the paraxial approximation, is equal 



to 



*(p,q) = Cexp(-^|p + q| 2 



(3) 
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where C is a normalization constant ensuring 
J dq J dp\^> (p, q) | 2 = 1. After integration along 
the z-axis one obtains 



*(p,q) = C 1 



-exp 
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(4) 



where erf refers to the error function. Notice that Eq. 
Q is similar to the one describing the joint spectrum of 
photon pairs in the frequency domain, when the spatial 
degree of freedom is omitted EH W 



Since all the configuration parameters that define the 
down conversion process show rotational symmetry along 
the propagation direction z, the joint amplitude given by 
Eq. (|4]) can be written as 



* (p, q) 



Bi (p, q) e 



H(<Pp—<Pq) 



(5) 
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Here, we have made use of polar coordinates in 
the transverse wave-vector domain for the signal, 
p = (p cos (p p ,p sin ip p ), and idler photons q = 
(q cos ip q , q sin ip q ) , where (p p<q are the corresponding az- 
imuthal angles, and p, q are the radial coordinates. The 
specific dependence of the Schmidt decomposition on the 
azimuthal variables tp p and ip p reflects the conservation 
of orbital angular momentum in this SPDC configuration 
21], so that a signal photon with OAM winding number 
+1 is always accompanied by a corresponding idler pho- 
ton with OAM winding number —I. The probability of 
such coincidence detection for each value of I is the spi- 
ral spectrum [22[ of the two-photon state, i.e., the set 
of values Pi = J pdpqdq\Bi(p,q)\ 2 . Recently, the spiral 
spectrum of some selected SPDC configuration have been 
measured t 23j. 

The Schmidt decomposition [3, of the spiral func- 
tion, i.e., Bi{p, q) = Y,n=a y/^nlfnl(p)gnl(q), is the tool to 
quantify the amount of entanglement present. A„; are the 
Schmidt coefficients (eigenvalues), and the modes /„; and 
g n i are the Schmidt modes (eigenvectors). Here we obtain 
the Schmidt decomposition by means of a singular-value 
decomposition method. Once the Schmidt coefficients 
are obtained, one can obtain the entropy of entangle- 



ment as E = — ^2 nl \ n i log 2 X n i- An estimation of the 
overall number of spatial modes generated is obtained 
via the Schmidt number K = 1/ J2 n i *nl > which can be 
interpreted as a measure of the effective dimensionality 
of the system. Finally, the spiral spectrum is obtained as 
p i = En ^u- 

For the sake of comparison, let us consider first the 
usual case of a QPM crystal with no chirp, i.e., a = 
0/xm -2 , and length L = 20 mm, pumped by a Gaussian 
beam with longitudinal wavenumber k p — 15.7 fim^ 1 and 
pump beam waist wq = 200 p,m. In this case, the inte- 
gration of Eq. ^ leads to (26J 



*(p.q) = Cexp 



Wo. 



Lip- 



4fc„ 



(6) 

The Schmidt coefficients are plotted in Fig. [SJa) , and the 
corresponding spiral spectrum is shown in Fig. [3](a) . The 
main contribution to the spiral spectrum comes from the 
spatial modes with I = 0. The entropy of entanglement 
for this case is E = 5.1 ebits and the Schmidt number is 
K = 18.3. 

Nonzero values of the chirp parameter a lead to an in- 
crease of number of generated modes, as it can be readily 
seen in Fig. ^h) for a — 10 x 10~ 6 fj,m~ 2 . This broaden- 
ing effect is also reflected in a corresponding broadening 
of the spiral spectrum, as shown in Fig. [3fb). Indeed, 
Fig. Hla) shows that the entropy of entanglement in- 
creases with increasingly larger values of the chirping pa- 
rameter, even though for a given value of Wq, its increase 
saturates for large values of a. For wq — 200 pm and 
a = 10 x 10~ 6 fim^ 2 , we reach a value of E = 13.47 
ebits. On the contrary, the Schmidt number K rises lin- 
early with a, as can be observed in Fig. [U^b), for all 
values of wq. For sufficiently large values of wq and a, 
K reaches values of several thousands of spatial modes, 
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FIG. 3. The spiral spectrum Pi for (a) a — /im~ and (b) 
a — 10 x 10~ 6 fim~ 2 . The nonlinear crystal length is L = 20 
mm and the pump beam waist is wo — 200 \xm. 



FIG. 4. (a) The entropy of entanglement E and (b) the 
Schmidt number K as a function of the chirp parameter a 
for wo — 40 /im (solid black line), wo = 100 /im (dashed blue 
line) and wo = 200 fim (dotted-and-dashed red line). 



i.e. K = 8097.76 for the same wo and a. For large val- 
ues of E, a further increase of E requires an even much 
larger increase of the number of spatial modes involved, 
which explain why an increase of the number of modes 
involves only produces a modest increase of the entropy 
of entanglement. 

For the sake of comparison, when considering fre- 
quency entanglement, the entropy of entanglement de- 
pends on the ratio between the bandwidth of the pump 
beam (typically B p ~ 5 MHz) and the bandwidth of the 
down-converted two-photon state (Bd c ) [27| . For type II 
SPDC, one has typically values of E ~ 1 - 2 Q- In- 
creasing the bandwidth of the two-photon state, one can 
reach values of B^c > 1000 THz, therefore allowing typi- 
cal ratios greater than B^j B p 10 s , with E > 25 [281 ] . 

In conclusion, we have presented a new way to increase 
significantly the amount of two-photon spatial entangle- 
ment generated in SPDC by means of the use of chirped 
quasi-phase-matching nonlinear crystals. This opens the 
door to the generation of high entanglement under var- 
ious experimental conditions, such as different crystal 
lengths and sizes of the pump beam. 

QPM engineering can also be an enabling tool to gener- 
ate truly massive spatial entanglement, with state of the 
art QPM technologies potentially allowing to reach 
entropies of entanglement of tens of ebits. Therefore, the 



promise of reaching extremely high degrees of entangle- 
ment, offered by the use of the spatial degree of freedom, 
can be fulfilled with the scheme put forward here. The 
experimental tools required are available nowadays. 

The shaping of QPM gratings are commonly used in 
the area of non-linear optics for multiple applications 
such as beam and pulse shaping, harmonic generation 
and all-optical processing [29j . In the realm of quan- 
tum optics, its uses are not so widespread, even though 
QPM engineering has been considered, and experimen- 
tally demonstrated, as a tool for spatial 3^, 31 1 and fre- 
quency [ij} control of entangled photons. In view of the 
results obtained here concerning the enhancement of the 
degree of spatial entanglement, it could be possible to de- 
vise new types of gratings that turn out to be beneficial 
for other applications in the area of quantum optics. 
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